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Abstract: We point out that perturbation theory in conjunction with the renormahzation 

group {RG) puts a severe constraint on the structure of the large-A^ non-perturbative 

CN \ glueball propagators in SU{N) pure YM, in QCD and in J\f = 1 SUSY QCD with massless 

Qv I quarks, or in any confining asymptotically- free gauge theory massless in perturbation theory. 

^£ , For the scalar and pseudoscalar glueball propagators in pure YM and QCD with massless 

\^ . quarks we check in detail the i?G-improved estimate to the order of the leading and next- 

'sj" . to-leading logarithms by means of a remarkable three-loop computation by Chetyrkin et al. 

^f^ . We investigate as to whether the aforementioned constraint is satisfied by any of the scalar 

or pseudoscalar glueball propagators computed in the framework of the AdS String/ large- 

N Gauge Theory correspondence and of a recent proposal based on a Topological Field 

kS ' Theory underlying the large- A^ limit of YM. We find that none of the proposals for the 

H ■ 
Cu ■ scalar or the pseudoscalar glueball propagators based on the AdS String/ large-A Gauge 

Theory correspondence satisfies the constraint, but the Topological Field Theory does. 
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1 Introduction and Conclusions 



1.1 Introduction 



In the last years several proposals for the non-perturbative glueball propagators of QCD- 
like confining asymptotically-free gauge theories have been advanced, based on the AdS 
String/ large- A^ Gauge Theory correspondence [1] and more recently on a Topological Field 
Theory (TFT) underlying the large- A^ limit of the pure Yang-Mills (YM) theory [2] [3] [4] 
[5] [6]. 

Above all these proposal aim to elucidate, at least in the large- A^ limit, the most fun- 
damental feature of the infrared of large- A^ QCD-like confining asymptotically-free gauge 
theories, i.e. the existence of a mass gap in the pure glue sector, as opposed to the massless 
spectrum of gluons in perturbation theory. 

However, these proposals predict a variety of spectra for large- A^ QCD different among 
themselves, asymptotically quadratic for large masses [1] [7] [8] [9] [10] or exactly linear 
[11] [3] [4] [5] ^ in the square of the glueball masses and in general do not agree about the 
qualitative and quantitative details of the low energy spectrum but for the existence of the 
mass gap. 

In view of the importance of the problem that these proposals aim to answer and in 
order to discriminate between the various proposals it is worth investigating if there is any 
constraint that we know by the fundamental principles of any confining asymptotically-free 
gauge theory that any supposed answer for the non-perturbative glueball propagators has 
to satisfy. 

In fact, we do know with certainty the implications of the asymptotic freedom for the 
large-momentum asymptotic behavior of any gauge invariant correlation function. 

In this paper we do not discuss at all the theoretical justification of the various proposal 
that we examine, leaving it to the original papers. We limit ourselves to check whether or 
not the constraint that follows by the asymptotic freedom and by the renormalization group 
in the ultraviolet (UV) is satisfied by any given proposal. Indeed, the importance of this 
constraint has been pointed out since the early days of large- A^ QCD [12], see also [13]. 

In fact the purpose of this paper is threefold. 

1.2 Implications of the renormalization group and of the asymptotic freedom 

Firstly, in sect. (2) we point out that perturbation theory in conjunction with the renormal- 
ization group (RG) severely constraints the asymptotic behavior of glueball propagators 
in pure SU{N) Yang-Mills, in QCD and in AA = 1 SUSY QCD with massless quarks, or 
in any confining asymptotically-free gauge theory massless to every order of perturbation 
theory. 



^Exact linearity in the TFT refers to the joint large-A'' spectrum of scalar and pseudoscalar glueballs. 
The TFT in its present formulation does not contain information about higher spin glueballs. 



Indeed, we show in this paper, on the basis of RG estimates, that the most fundamental 
object involved in the problem of the mass gap ^, the scalar [S) glueball propagator in any 
(confining) asymptotically-free gauge theory with no perturbative physical mass scale, up to 
unphysical contact terms, i.e. distributions supported at coinciding points, has the following 
universal, i.e. renormalization-scheme independent, large-momentum asymptotic behavior: 
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Analogously for the pseudoscalar (PS) propagator: 
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and for a certain linear combination as well, the anti-seldual (ASD) propagator: 
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terms. The coincidence of the asymp- 



totic behavior, up to the overall normalization constants that are computed in sect. (3), 
Cs,Cps,Casd, is due to the coincidence of the naive dimension in energy, 4, and of the 
one-loop anomalous dimension, 7(51) = —2/305^ + ' ' ' j of these operators deprived of the 
factors of ^^ or of g^. Euclidean signature is always understood in this paper unless 
otherwise specified. 



^The lightest glueball is believed to be a scalar in pure YM and in the 't Hooft large- A^ limit of QCD. 



1.3 Perturbative check of the RG estimates 

Secondly, in sect. (3) we check the correctness of our RG estimate on the basis of an exphcit 
very remarkable three-loop computation performed by Chetyrkin et al. [14] [15] in pure 
SU{N) YM and in S'C/(3) QCD with nj massless Dirac fermions in the fundamental 
representation. For example, we show that in pure SU{N) YM Chetyrkin et al. result [14] 
[15] can be rewritten by elementary methods as: 
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where a and a are two scheme-dependent constants that are defined in sect. (3.5) and 

is the 't Hooft coupling constant in the MS scheme. In Eq.(1.4) the terms that depend on 

2 

g{-r7 — ) correspond in the coordinate representation to distributions supported at coincident 

Ws 
points (contact terms), and therefore they have no physical meaning. Remarkably, the 

correlator without the contact terms does not in fact depend on the arbitrary scale // 

(within 0{g^) accuracy) as it should be. The running coupling constant 5xrc( a^ — ) occurs 

in Eq.(1.4) with two-loop accuracy and it is given by: 
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Therefore, the perturbative computation furnishes an expansion of the correlator in powers 
of g^-—(^fj,) and of logarithms. This expansion has been rearranged by elementary methods 
in terms of the two- loop running coupling g^^-f-^{-fM — ) in Eq.(1.4). 

At this point our basic strategy to check the RG estimates of sect. (2) consists in 
substituting in Eq(1.4) instead of Eq.(1.5) the i?G-improved expression for g^—^(-j^ — ) given 
by: 
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The MS scheme is indeed defined [16] in such a way to cancel the term of order of -^ — - — : 
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that would occur in Eq.(1.6) in other schemes. By subtracting the unphysical contact 
terms and by substituting the i?G-improved two- loop asymptotic expression for g^—-[j^- 
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it follows the actual large-momentum scheme-independent asymptotic behavior of Eq.(1.4): 
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as opposed to the perturbative behavior that would follow by Eq.(1.5). The asymptotic 
result in the other cases is checked similarly. 

1.4 AdS/harge-N Gauge Theory correspondence and disagreement with the 
RG estimates 

Thirdly, in this subsection and in the next one, we inquire whether the large- A'^ non- 
perturbative scalar or pseudoscalar propagators actually computed in the literature agree 
or disagree with the RG estimate. 

We find, to the best of our knowledge, that all the scalar propagators presently com- 
puted in the literature in the framework of the AdS String/ large- A^ Gauge Theory corre- 
spondence disagree with the universal asymptotic behavior. 

We should mention that the comparison of the asymptotics of the scalar glueball prop- 
agators in the AdS approach with YM or with QGD at the lowest non-trivial order of 
perturbation theory has been already performed in [17] [18] [19] [20], but with somehow 
different conclusions. The reasons is that in [17] [18] [19] [20] the comparison has been 
performed only with the one-loop result for the scalar glueball propagator, i.e. only with 
the first term in Eq.(3.8), that is conformal in the UV. No higher order of perturbation 
theory and no RG improvement has been taken into account in the comparison, as instead 
we do in this paper. 

Here we enumerate the models based on the AdS /Gauge Theory correspondence for 
which we could find explicit computations of the scalar glueball propagator in the literature. 

In the Hard Wall model (Polchinski-Strassler background [9] in the so called bottom-up 
approach) : 
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where Ki,Ii are the modified Bessel functions [17]. The asymptotic behavior [17] is con- 
formal in the UV: 
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with p = \/p^. Indeed, as recalled in appendix A, in the coordinate representation: 
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and, as observed in [17] [18] [19] [20], it matches the one-loop large- A^ QCD result for the 
perturbative glueball propagator displayed in the first term of Eq.(3.8). Nevertheless, it 
disagrees by a factor of (logp)" with the correct asymptotic behavior in Eq.(l.l). 

The Soft Wall model (bottom-up approach) [11] implies the same leading conformal 
asymptotic behavior [17] [18] [19] ]20] in the UV for the scalar glueball propagator: 



/ {trF'^{x)trF^{0)) e-^P'^d^x ~ -/ 
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that therefore disagrees in the UV by the same factor of (logp)^. 

A more interesting example of the AdS string / large- A^ Gauge Theory correspondence 
from the point of view of first principles applies to the cascading M = 1 SUSY YM 
theory (top-down approach) [21] [22], because in this case the correct asymptotically-free /3 
function of the cascading theory is exactly reproduced in the supergravity approximation 
in the Klebanov-Strassler background [21] [22]. Nevertheless, the asymptotic behavior of 
the scalar correlator is [23] [24]: 
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that disagrees by a factor of (logp)^ with the correct asymptotic behavior in Eq.(l.l). 

1.5 Topological Field Theory and agreement v^ith the RG estimates 

Finally, in sect. (4) we prove that in the large- A^ limit of pure SU{N) YM the ASD glueball 
propagator computed in [5] [4] [3] ^: 
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al3 al3 k=l ^ W 



agrees with the universal RG estimate in Eq.(1.7). 

Since the proposal for the TFT underlying large- A^ YM is recent and not widely known 
we add here a few explanations, but for the purposes of this paper the reader can consider 
Eq(1.8) just as a phenomenological model factorizing the ASD glueball propagator on a 
spectrum linear in the masses squared with certain residues. 

Yet, to say it in a nutshell, the rationale behind Eq(1.8) is as follows. In [2] [4] [6] it is 
shown that there is a TFT trivial [4] [6] at A^ = oo underlying the large- A^ limit of YM. 



^We use here a manifestly covariant notation as opposed to tiie one of tiie TFT employed in [5] [4] [3]. 



At A^ = CX3 the TFT is localized on critical points [3] [4]. However, at the first non-trivial 

jf order the ASD propagator of the TFT arises computing non-trivial fluctuations around 

the critical points of the TFT [3] [5]. 

In Eq.(1.8) F^a is the anti-selfdual part of the curvature F^p = daAp — dpA^ + 

i-j=[Aa,Ajj\ with the canonical normalization defined in Eq.(2.8), K-^ is the renormal- 

ization group invariant scale in the scheme in which it coincides with the mass gap and 

9k = di^T- = k) is the 't Hooft running coupling constant at the scale of the pole (in 

w 
Minkowski space-time) in the scheme defined in [2], that is recalled in sect. (4). In fact, 

the analysis of the UV behavior of Eq.(1.8) has already been performed at the order of 

the leading logarithm occurring in Eq.(1.3) in [5]. Here we go one step further comparing 

Eq.(1.8) with Eq.(1.3) at the order of the next-to-leading logarithm. Our basic strategy 

to obtain the large momentum asymptotics of Eq.(1.8) is as follows. We write the RHS of 

Eq.(1.8) as a sum of physical terms and contact terms according to [5]: 

k=l ^ W k=l ^ W fc=l 

The first sum contains the physical terms that in Minkowski space-time carry the pole sin- 
gularities, while the second sum contains the contact terms, that we ignore in the following. 
We now consider only the physical terms and to find the leading UV behavior we use the 
Euler-McLaurin formula according to the technique first introduced by Migdal [25] ^ and 
employed in [5] : 

E Gkip) = / G,{p)dk -Y,-i di-'G,{p) (1.10) 
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where Bj are the Bernoulli numbers. In our case the terms proportional to the Bernoulli 
numbers involve negative powers of p and they are therefore suppressed with respect to the 
first term which behaves as the inverse of a logarithm, so that we ignore them as well. We 
have: 
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where we have used the i?G-improved asymptotic behavior for large k of the running cou- 



We understand that Migdal technique has been known to him for decades. 



pling constant (7^ at the scale of the k-th pole, i.e. on shell (in Minkowski space-time): 
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The constant c is related to the scheme that occurs in the non-perturbative calculation 
[2] [3] [4] [5] . The actual value of c is not relevant in this paper since we study only the 
universal asymptotic behavior. In sect. (4) we compute the universal leading and next-to- 
leading behavior of the integral in Eq.(l.ll) and the result is: 
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1.6 Conclusions 



The preceding result, for the ASD glueball propagator computed in the TFT underlying 
large-A^ pure YM, agrees perfectly in the large-A^ limit with the universal part of the 
renormalization group improved expression of the perturbative result Eq.(1.7). 

The agreement is due to the conspiracy between the residues of the poles, that are 
proportional to the fourth power of the coupling constant renormalized on shell times the 
fourth power of the glueball mass at the pole, and the exact linearity of the joint scalar and 
pseudoscalar spectrum of the square of the mass of the glueballs in the ^45-0 correlator of 
the TFT. 

To the best of our knowledge this is the only non-perturbative result for the scalar or 
pseudoscalar glueball propagator proposed in the literature that agrees with large- A^ YM 
perturbation theory and the renormalization group. 

While this agreement is not by itself a guarantee of correctness of Eq.(1.8) it deserves 
further investigations, both at theoretical level and of further checks. 

Besides, our analysis shows that the AdS/harge-N Gauge Theory correspondence in 
any of its present strong coupling incarnations, the bottom-up or the top-down approach, 
for which scalar glueball propagators are available in the literature, does not capture, not 
even approximatively, the fundamental ultraviolet feature of YM or of QCD or of any 
large- A^ confining asymptotically-free gauge theory in the pure glue sector. 

While this conclusion is certainly known to some experts (see for just one example [22]), 
we think that it is not widely recognized that constructing theories that are conformal in 



the ultraviolet, as the Hard or the Soft Wall models, or even with the correct beta function 
but in the strong coupling phase, as the Klebanov-Strassler supergravity background, is 
not at all a good approximation for the correct result in the ultraviolet. In this paper, 
for the first time with leading and next-to-leading logarithmic accuracy, we have computed 
quantitatively the measure of the disagreement. 

Finally, given the disagreement between the propagators of the TFT and the propa- 
gators of the AdS /Ijaxge-N Gauge Theory correspondence in the infrared for the first few 
lower- mass glueballs, a careful critical analysis of the two approaches at level of numerical 
lattice data is needed, and also at theoretical level of further constraints arising by the 
OPE and by the low-energy theorems of Shifman-Vainshtein-Zakharov [SV Z). 

2 Renormalization group estimates on the universal behavior of corre- 
lators 

2.1 Definitions 

The SU{N) pure YM theory is defined by the partition function: 

■ JtrF'\x)d'^x 
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where we use the simplified notation trF'^{x) = X^^fl ^^(-^ofl)- Introducing the 't Hooft 
coupling constant g [26]: 

g^ = g^MN (2.2) 

the partition function reads: 

Z= (vAe-^^'-''^^^''^ (2.3) 



According to 't Hooft [26] the large- A^ limit is defined with g fixed when N — )■ oo. 

For the structure of large- A/^ glueball propagators see [12] and for reviews of the large- A 
limit see [13] and [27]. The normalization of the action in Eq.(2.1) corresponds to choosing 
the gauge field Aa in the fundamental representation of the Lie algebra, with generators 
normalized as: 

tr {ft^) = h"^ (2.4) 

In Eq.(2.1) Fajj is defined by: 

F„/3(x) = do^Ap - dpAo, + i[Ao,,Ap\ (2.5) 

We refer to the normalization of the action in Eq.(2.1) as the Wilsonian normalization. 
Perturbation theory is formulated with the canonical normalization, obtained rescaling the 



field Aa in Eq.(2.1) by the coupling constant gvM = -f^' 

Aa ^ <7yM^r (2.6) 

in such a way that in the action the kinetic term becomes independent on g: 

]- f tr{F\A''^'')){x)d^x (2.7) 

where: 

Fai3 = df.A'^r - d^Af"^ + igvM [AT, Af^] (2.8) 

From now on we will simply write F^p for the curvature as a function of the canonical field, 
without displaying the superscript can. 

2.2 A short summary of perturbation theory and of the renormaUzation group 

We recall the relation between bare and renormalized two-point connected correlators of 
a multiplicatively renormalizable gauge-invariant scalar operator O of naive dimension in 
energy D: 

G(^\p,^^,g{^^)) = Zl{-,g{A))G^^\p,A,g{A)) (2.9) 

(2) 

where Gq is the bare connected correlator in momentum space computed in some regu- 
larization scheme with cutoff A and /i is the renormalization scale: 

Gj')(p,A,5(A)) = y"(O(x)O(0))_,„e^^-d^x^(O(p)O(-p))_„ (2.10) 

Since YM or QCD with massless quarks or A/" = 1 SUSY YM with massless quarks is 
massless to every order of perturbation theory and since O has naive dimension D we can 
write: 

G^'\p,f,,gif,))=p'^-'G^^l{^,gi^t)) (2.11) 

(2) 

where G]^£ is a dimensionless function. The Callan-Symanzik equation for the dimension- 
less two-point renormalized correlator expresses the independence of the bare two-point 
correlator with respect to the subtraction point fi: 
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where we have defined the beta function with respect to the renormahzed couphng g{fi): 

m = -7^ , ,,, (2.14) 

and the anomalous dimension: 

loia) = — ?n (2.15) 

91og^ A,a(A) 

Eq.(2.13) can be rewritten taking into account the dependence of G^^ on the momentum 
P = yP^'- 



d 



^ d log p 
The general solution of Eq.(2.16) is: 



/3(5)^-27o(^))G'g[(^,5(/.)) = 
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/i // /x 

The running coupling g{-,g{fi)), that we briefly denote by g{p), solves: 
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with the initial condition g{l,g{ii)) = g{n). 

The multiplicative renormahzed factor Zo{^^g{^)) satisfies: 
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and from now on it is thought as a finite dimensionless function of g(fj,) and g{p) only 
Zo{9{p),9{lJ'))- 

rg(p) 7o(g) j„ 

Zo = e-'9('^) P^a) "3 (2.20) 

Eq.(2.17) expresses the solution of the RG equation as a product of a RG invariant (RGI) 
function Qq of g{p) only and of a multiplicative factor Z^ that is determined by the anoma- 
lous dimension "foid) ^•iid by the beta function P{g). Go a-i^d Zq can be computed order 
by order in renormahzed perturbation theory. 

From Eq.(2.18), that represents the coupling constant fiow as a function of the mo- 
mentum, we obtain the well-known behavior of the iJG-improved 't Hooft running coupling 
constant with one- and two-loop accuracy, starting from the one- and two-loop perturbative 
beta function: 

m = -Pog^ - fiig^ + ■■■ (2.21) 
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where for pure YM: 



A " ' 
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A = ^77^ (2-22) 



With two-loop accuracy we get 
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This is the well known actual UV asymptotic behavior of the running coupling constant. 
The function Qo in Eq.(2.17) is not known from general principles but can be computed 
in perturbation theory as a function of g{fj,) and then expressed in terms of g{p), since Qo 
is RGI. Similarly, we can evaluate Zq using again the one-loop or two-loop perturbative 
expressions for l3{g) and 70(5): 

lo{g) = -70(0)5^ - 7i(o)9^ + • • • (2-24) 

With one-loop accuracy: 
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and with two-loop accuracy we have: 
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(7o(ci)7i(ci)/3o/5i - (7o(ci)/5i 



log log ^ 



^m<^gif 



(2.26) 



In evaluating the last two expressions we have used the two-loop i?G-improved expression 
for g[p) given by Eq.(2.23). From the two- loop -RG-improved expression in Eq.(2.26) it 
follows that the leading and next-to-leading logarithms for Zq are determined only by /3o, 
/3i and by 70(0), that are in fact universal, i.e. scheme independent. Indeed, the two-loop 
coefficient of the anomalous dimension 71(0) does not occur in the ffist log log - term, but 
only in terms that have a subleading behavior as powers of logarithms. Keeping only up 
to the next-to-leading term in Z"^, we obtain for the universal logarithmic behavior of the 
dimensionless two-point correlator: 






^2/3olog^ 



/3i log log ^ ' 
'" 2/30^ log f 



T0(O) 
/30 



Qo{g{p)) 



(2.27) 



Thus our aim, in order to get asymptotic estimates, is to determine the one-loop coefficient 
of the anomalous dimension 70(0) and the RGI function Qq for our operators O. 

2.3 Anomalous dimension of trF"^ and of trFF 

The operator B^trF"^ is proportional to the conformal anomaly, that is the functional 

derivative with respect to a conformal rescaling of the metric of the renormalized effective 

action that must be RGL Therefore ^-^trF"^ is RGI as well. Hence its anomalous di- 

9 

mension vanishes and, using the notation of the previous section, the form of its correlator, 
ignoring possible contact terms that will be taken into account in sect. (3), is: 



9 



' P^Qpis)p2 



{9{P)) 



(2.28) 
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On the other hand trF'^ is not RGI and therefore its correlator is: 

Gpl (p, /i, gifj-)) = p^Qf-^ igip))Zl2 (77, ^(m)) 



Ai 



(2.29) 



(2) (2) 

Since the relation between G^(g) ^{p^ fij g{fi)) and Gj^2{Pf f^^g{f^)) is: 



G^g,^^{p,H,g{f,))=(^) Gg](p,M,5(M)) 

9 \ y / 



(2.30) 



it follows that ( ^-^ J and Z'^{-,g{fi)) must combine in such a way to obtain a function of 
ff(p) only: 



?2 ,„„ „,,^i^;+ ^vi„„ .,,v„;-r, +u^ property of ^-i^ ■'■'»■ ^^ 
being RGI. Its two-point correlator must indeed satisfy the equation: 



To find the anomalous dimension ^p2 of trF'^ we exploit once again the property of ^-^trF'^ 



d 



d 



dlogp 



.(2) 



f^(9)jr-^ G)Z,(p,f,,g{^,))=0 



F2~ 



(2.32) 



where the last term occurs because we are considering the complete correlator and not its 
dimensionless part. Using Eq.(2.30) we find the anomalous dimension of trF"^: 



d 



dlogp 



m 



d 



/5(5)V^(2) 



GUip,fJ-,gifi)) 







g J dlogp \ g J dg 



d 



d 



dlogp dg dg \ g J 

From this equation it follows: 

d fl3{g) 



Gf,{p,fi,g{fi)) = 
Gp,{p,fi,g{fi)) = 



-ip2{g)=g- 

dg \ g 

With two-loop accuracy this expression reads: 

7irF2(ff) = -2/3o5'-4/3i/ + 



(2.33) 



(2.34) 



(2.35) 



14 



72 (P 
Jill J' Lilt; iiiDL Lc:iiii, vv^ \^a,LL viciivc; Liic ca.^1 coDiuii L\JL z. — 

accuracy: 



Keeping only the first term, we can derive the expression for Z'^(-,g(fj:)) with one-loop 



Z\^,9{f^))-^, (2.36) 



/3(9), 



Finally, the correlator of ^-^trF , with one-loop accuracy, is: 

G'-g,^^{p,fi,g{f,)) =p'Pig\^^)p^^gM9{p)) =p'Pi9Hp)gF<9{p)) (2-37) 

We can repeat similar calculations for the operator trFF in order to compute its anomalous 
dimension, using the property of g^trFF being RGL Indeed g^trFF is the density of the 
second Chern class or topological charge. The Callan-Symanzik equation is: 



dlogp 



m^g-^)G%piP.^^,9{^^)) 



' /3(.)#-4 



g'G%{p,f,,g{f,)) 



d log p dg 

from which we obtain the anomalous dimension of trFF: 



(2.38) 



7^^(5) = 2^ = -2/305^ - 2/3i/ + • • • (2.39) 

9 

We notice that while the one-loop anomalous dimensions of trF'^ and of trFF coincide, the 
two- loop anomalous dimensions are different. This means that the operator trF~ has a well 
defined anomalous dimension only at one loop, in agreement with the fact that it belongs to 
the large- A^ one-loop integrable sector of Ferretti-Heise-Zarembo [28]. Therefore, only the 
universal part of its correlator, that is determined by the one-loop anomalous dimension 
and by the two-loop /3 function, can be meaningfully compared with the non-perturbative 
computation in Eq.(4.1). 

2.4 Universal behavior of correlators 

Knowing the naive dimension D and the anomalous dimension of a (scalar) operator Od, 
we can write the asymptotic form for p >> ^ of its correlator obtained by the RG theory. 
Indeed, as we recalled in sect. (2.2), assuming multiplicative renormalizability, the RG- 
improved form of the Fourier transform of the correlator is given by: 

G(2)(p2)= f{On{x)ODm,,^^e^P-d'x = p'^-'ga,{9{p))ZljP,g{f^)) (2.40) 

J /^ 

where the power of p is implied by dimensional analysis, Qod i^ ^ dimensionless function 
that depends only on the running coupling g{p), and Zq is the contribution from the 
anomalous dimension. But in fact in general the correlator of Od is not even multiplicatively 



renormalizable because of the presence of contact terms. These terms would affect the UV 
asymptotic behavior, but they are non-physical and therefore they must be subtracted. 
In fact, they spoil the positivity of the correlator in Euclidean space in the momentum 
representation, that is required by the Kallen-Lehmann representation (see the comment 
below Eq.(3.21)). 

In the coordinate representation of the correlator, for x ^ 0, contact terms do not occur. 
Therefore, a strategy to avoid that contact terms interfere with the RG improvement is to 
pass to the coordinate scheme [29], where the correlator is multiplicatively renormalizable, 
to compute its i?G-improved expression, to go back to the momentum representation, and 
eventually to subtract the contact terms. 

In the coordinate representation for x ^ the solution of the Callan-Symanzik equation 
reads: 

G§l{x) = (Oz5(rE)O^(0)),„„„ = {^fgoMx))Zl,Jx^^,g{^,)) (2.41) 

with x = va?, where we have denoted by g[x) the running coupling in the coordinate 
scheme [29] and by an abuse of notation we have used the same names Q and Z for the 
RGI function and renormalization factor in the coordinate and momentum representation. 
The function Q{g{p)) can be guessed at one loop since the correlator must be confor- 
mal at one loop in perturbation theory, that implies Qii-ioop){9{^)) = const. Indeed, for 
composite gauge invariant operators the renormalization of the coupling occurs starting at 

two loops only. Hence: 

P P 

G(i-ioop){-) = const log - (2.42) 

Indeed in appendix A we show that Jp'^logj^e^^'^d'^p = const{^x) that is conformal 
in the coordinate representation. The explicit dependence on ^u, that contradicts RG in- 
variance in the momentum representation, is due to the fact that the correlator in the 
momentum representation, as opposed to the coordinate representation, is not really mul- 
tiplicatively renormalizable because (scale dependent) contact terms arise. This can be 
understood observing that in the coordinate representation for x 7^ the one-loop correla- 
tor is independent on the scale // but it is not an integrable function, in such a way that its 
Fourier transform needs a regularization, that introduces the arbitrary scale ^u. 
Naively, we can already derive the leading UV asymptotic behavior: 

Ggi(p^)~p--log^(fjM) "' -p'-~\9\p)r^-' (2.43) 

where we have used the fact that g^ij)) ~ t — tpt- It easy to check that for D = A and 
1o(Od) ~ 2/^0 this estimate coincides with Eq.(l.l)-Eq(1.3). 

However, this estimate assumes multiplicative renormaliz ability in the momentum rep- 
resentation and it does not take into account the occurrence of contact terms in the mo- 
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mentum representation of the correlators. 

Nevertheless, in the next section we confirm by direct computation that after subtract- 
ing the contact terms the actual behavior of the scalar and of the pseudoscalar correlator 
agrees with the estimate in Eq.(2.43). 

3 Perturbative check of the universal behavior of correlators 

In this section we obtain the explicit form of the three-loop correlators of trF"^ and of 



trFF starting from their imaginary parts that have been computed in [14] [15] in the MS 
scheme. The MS scheme can be defined as the scheme in which the two-loop i?G-improved 



running coupling does not contain 



l°g'fe 



contributions [16]. More precisely, we consider 



the equation for the running coupling constant that follows from the two-loop /3 function: 

/■^(P) _dg_ _ 1_ 



1 P 



+ ^log{g{p))+C + --- 



(3.1) 



where C is an arbitrary integration constant and A^ is the RGI scale in a generic scheme 
s. The value of C in the MS scheme is chosen in such a way to cancel the — - — 3- term in 



the solution: 



log^ 



/3o log ^ 


1- 


/3ilog(/5ologfj) ^ c 
PI log^ ' logf^ 


C=§log(/3o) 
Po 







(3.2) 



The result reported in [14], for trF"^ in the SU{3) YM theory, is: 



lm{trF\p)trF\-p))^^^^ = ^/|l + 

vr 



vr 
37631 
96 
363 



21 



73 11 , p 
^C(2)-?C(3) 



2817 , p^ 
——- log — + —T 
16 n'^ 16 



1 2P 
log -^ 



(3.3) 



where a^ = -^^ and aj^ is a^ in the MS scheme. Firstly, we want to find from Eq.(3.3) 
the result for the SU{N) YM theory and we want to express the result in terms of the 
't Hooft coupling in the MS scheme ffjjg- In fact, this operation is quite easy since it is 
known, and it can be checked in [15], that at this order of perturbation theory the rank of 
the gauge group enters the result only through the Casimir factor Ca = N. Therefore, to 
obtain the general result it is sufficient to divide by 3 and to multiply by A^ the coefficient 
of ctf^ and to divide by 9 and to multiply by A^^ the coefficient of djjq ■ The factors of A^ 
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and of N"^ are then absorbed in the definition of the 't Hooft coupling constant. We obtain: 



lui{trF\p)trF\-p))^^^„ 



47r 



73 



2 ^,. .. log 



3(47r)2 3(47r)2 ^^ ^'^ 



+ 9^(lA 



37631 



242 



110 



54(4^)4 3(4^)4^^^^ (47r)4^^^^ 



313 , p2 121 , oP 

log — + TTTTZTJ log -^ 



(4^)4 ^^2 3(4^)4 



(3.4) 



From Eq.(3.4) we derive the complete expression of the correlator, assuming the correlator 
in the form: 



{trF'^{p)trF^{-p))^ 



N^-l 4i P 



1 + 5|7s(/^) ( /o - /3o log ^^2 



El 



+ <?1^(/^) ( /l + /2 log -2 + /3 log^ -2 



(3.5) 



We extract the imaginary part of Eq.(3.5) that arises from the imaginary part of the loga- 

2 2 

rithm in Minkowski signature, log(— ■^) = log ^ — in. We obtain: 



Im {trF'^{p)trF^{-p)) 
{N^ - 1) 



A-K 



-P 



conn 
,2/,.\ , r f f „2\„4/ 



l + /05l/^) + (/i-/3vr^)5'(M) 



2 2 2 " 

2/3o5'(/u) log ^ + 2/2/(/i) log ^ + 3/3/(m) log' ^ 

;U^ /X^ /i^ 



(3.6) 



Finally, comparing Eq.(3.4) with Eq.(3.6) we determine the values of the coefficients /j: 

73 



/o 



3(47r)2 
2 _ 37631 242 
~ ^^4 3~ 



-2/3o 

2/2 
3/3 



11 



3(47r)2 
313 



/2 



C(2)-110C(3))^ 



313 



(47r)4 

121 

4^/3 



3(47r 



2(47r)4 



9(47r) 



/. = (5^_noc(3)) ' 



54 



(47r)^ 



(3.7) 
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Therefore, the correlator is: 
{trF\p)trF\-p))^^ 



{N^ - 1] 
47r2 



■P log 



p^ 



l + 5'(/^) /o-/3olog 



/"^ 



+ /(/^)(/i + /2log^ + /3o'log2^ 



(3.8) 



Similarly, the imaginary part of the correlator oftrFF, already written in [15] for the gauge 
group SU{N), is: 



Im {trFF{p)trFF{-p))^ 



47r ^ \ TT [ V12 6 ^ iJ,^ 



+ ( 



«Ms(^) • 



A..r^_^C(2)-fc(3) 



vr L V 864 

1135 , p2 121 , 2 P^ 



(3.9) 



We obtain: 

{trFF{p)trFF{-p))^ 



(A^2 _ 1) 
47r2 



p log 



A.2 



1 + S'ItqI/") /o - ^0 log 






+ 5i75 (^) ( /i + /2 log ^ + /3o' lor' ^' 



/^" 



/^ 



(3.10) 



where: 



97 



3(4^) 



A = ,5]|^_noc(3)) > 



54 



(4vr) 



-2/3o 

2/2 



11 



3(4^)2 

1135 ~ 

4^/2 



3(47r 



1135 

'6(47r)4 



(3.11) 



3.1 Correlator of ^trF'^ in 5[/(7V) FM (two loops) 



gAT 



We now determine the UV asymptotic behavior for the correlators by employing their 
i?G-improved expression. Firstly, we recall that in every generic scheme labelled by a the 
relation between the coupling constant at two different scales is, with one-loop accuracy: 



1 1 ^^^^^P^ 



(3.12) 
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This relation is necessary to express the correlators in their i?G-improved form. As a 
starting simplified example we consider the two- loop expression of the correlator of jioj^trF'^ 
(for the moment we skip overall positive numerical factors in the normalization of the 
correlator): 



iv iV c 



2„4 4 



-/3oVff^ 



/^ 



MS"- l\2 

MS 



1 P' 



^* 



MS 



Po log -^ 



(3.13) 



This expression is renormalization group invariant with one-loop accuracy, since the factor 
^Pie2-j2 -g ^2^4 •£ ^g employ the one-loop /3 fund 
absorbed in a change of scheme. Indeed, defining: 



^Pie2-j2 -g ^2^4 •£ ^g employ the one-loop /3 function. The finite term fQg^-—[fi) can be 



5a (m) = 5|7c(^)(l + aglf^{^i)) 



MS"- 



(3.14) 



it follows: 



5a (/^) 



g|^(l + 2ag^ 



{^)+a'g^{fi))+0{g'') 



glrs^lA =gl{^i){l- agj^ifi) + a^g^i^^)) + 0(/) 

= gl{fi){l - aglif,) + 2a2<7a(/^)) + 0(/) 
5|^(/.) = 5.'(/z)(l - 2ag'M + 5a25,^(^)) + 0{g'') 



(3.15) 



We obtain for the correlator: 

9 2 

(/3o^trF2(p)/3o^trF2(-p)) 



-/3oVlog§ff,^(A.) 



2 2 " 

1 + (/o - 2a)<7^(/x) - Poglif^) log ^ + 0(5' log ^ 



To cancel the finite term it is sufficient to put: 



2 



(3.16) 



(3.17) 



Hence we obtain: 



{Po%trF\p)f3o'-trF\-p)) 

iV iV c 



P At 



-PoP log— 5a (m) 



l-f3oglil^)log 



1 2 ' 

\+0{g^\og^^ 



(3.18) 
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At this order the term in square brackets is precisely the renormahzation factor necessary 
to renormahze two powers of gad-i-)- We obtain: 

J V J V conn 

2 2 

- -/3ipSl{^i)9lip) log ^ (1 + 0(/ log ^)) (3.19) 

From Eq.(3.12) we express the logarithm in terms of the coupling constant: 

The correlator becomes: 

{Po^trF^{p)f3o^trF^{-p)) 

i V i V conn 

= Ap*(92(p) - 92(m)) (1 + Oig' log 4)) (3.21) 

The second term in the last line is in fact a contact term that has no physical meaning, 
therefore it may depend on the arbitrary scale /i since it must be subtracted anyway. The 
physical term is positive, despite the correlator that we started with was negative. This is 
an important feature, since a negative physical term would have been in contrast with the 
Kallen-Lehmann representation, that requires a positive spectral function. 



3.2 Correlator of ^trF"^ in SU{N) YM (three loops) 

We now consider the three-loop result Eq.(3.8), this time including also the correct normal- 
ization factors: 

ArF\pf^trF\-p)) 

iv iV conn 



1 N 4 / ^ P^ 1 P^ 



( P^ 

1 + djjs'yl^) ( /o - /5o log — 



+ 5fc(^)(/i + /2log§ + /3o2log2^ 



(3.22) 
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This correlator is not supposed to be RGI, because the factor of (HM)2 _ g4^-^_^ gd'^) is 
missing. We can ehminate the finite terms in the correlator by a redefinition of the scheme: 



albilJ') = 5|j5(/i)(l + ag^ifJ-) + bg^ifi)) 



Substituting we obtain: 



10\ 



(3.23) 



£trF\pf^trF\-p)) 

iV iV c 



-I 4 9 

I T) T) 

'^^ ~ l^)9lb{ij) (1 - '^agli^ifi) + (26 + 5a^)g^f,{n)) j^ log -^ 



47r2 

2 



P 



1 + /05'ab(^)(l - agabifJ-)) - /3ofi'a6(M)(l " agabif^)) log — + figabifJ-) 

+ f2gtbifJ') log -2 + l^htbif^) log^ -2 
/i /i _ 



.l-^)5afe(A^)^log^ 



1 + (/o - 2a)5i(/x) - /3ogl,{ti) log ^ + (/i + Sa^ + 26 - a/o)5^b(/i) 

A* 

2 2 " 

+ (/2 + 3/3oa)5^,(/i) log ^ + fSlglii,) log^ ^ 



/^ 



/^ 



(3.24) 



We eliminate the two finite terms choosing: 

/o 



/i + 5(^)2 + 26-^ = 



^b = -fl 



h 



(3.25) 



With this choice of a the coefficient of the g^ log ^ term becomes: 



/2 + 3/3oa = /2 + -/o/3o 



68 



3(4^) 



-2/3i 



(3.26) 



Therefore, the correlator reads: 



{UrF\p)UTF\-p)) 



P 



P 



'N 



N 



^i-iv^)^'^('^)4;^i°V 



1 - f^Ogabif^) log ^ - 2/3ig^fe(/i) log ^ + /3o5afe(A*) loL , 
/i /i /x 



(3.27) 
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We notice that the expression in square brackets is the two-loop Z factor determined by the 
anomalous dimension of trF"^ according to Eq.(2.35). The coefficient 1fi\ should become 
fi\ if we multiply the correlator in Eq.(3.27) by the factor of (l + g-gatil^)) i i'^ order to 
make the correlator RGI: 



c- 



Ngab 



trF\p)^^trF\-p)) 



Ngab 



1 



['^-^)pi9tb{iJ)-^{^+^jib{^j)y^og^ 



p-. 



/3o- 



/z2 



T) T) T) 

1 - /309afe(/^) log -^ - 2/3iC/^fe(/i) log -2 + /3o5a6(/^) log^ -2 
/X fl fl 

(l + |,2,(,))(l+ft,i(p))log^ 

T) T) T) 

1 - Poglbif^) log — - 2/3ig^fe(/i) log — + /SigihifJ') log^ — 

/Lt^ //^ /i^ 



(3.28) 



where we have multiplied and divided by (l + w-gabil^)) 1^ order to exploit the two-loop 



relation: 



Mi + ^glbip^og^ 



Po' 



1^' 



We now evaluate separately: 



9ab(.P) dabif^) 



(3.29) 



(1 + §;9Up)) 



00 ■ 



/3o 



^5afe(/^)+/5l5ab(/^)log— + 



^^^.(/^))(1 ^^-^-- «-^— ^' '^''^- 



/3o 



/^ 



/3i 







5^,(/i))+0(/log^) 



2 2 

l + /3iff^,(//)log^ + 0(<7^1og^; 



/x- 



(3.30) 
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Putting all together we get: 



(^^kAtrF\p)^^^trF\-p)) 



Ngab 



Ngab 

ab 
o2 



;i-^)/3o5a6(/^)^ 



4 
47^ 



;i + P,gl,{^^) log ^) (1 + ^gl,{^^)) (^ 



1 



1-/30506(^)10? 



/^ 



5ab(:P) 5a6(M) 

2p^gl,{^l)\og^ + PlgU^,)\o. 



) 



/i- 



1 ^ . . . . P^ 



/?: 



(l-T72)/305.^(/i)£2(l + ^5a.(/^))( 



1 



1 



iV2 



47r2 



/3o- 



ffafe(:P) 5afe(M) 
2 2 2 

T) 77 77 

1 - /3o5a6(^) log — - /3l5afe(^) log — + /3o5afe(/^) log^ — 
/i^ /i^ /X^ 



(3.31) 



The factor in square brackets in the last line is now precisely the renormalization factor for 
two powers of gat- Hence the correlator reads: 



^^iE!^trF\p)^^trF'' 



Ngab 



1 



Ngab 

P Ji 



-P)) 



/3: 



^ ~ i;H)(^0JZ23ab{lA9ab{p){^ + -^9ab{^J)) 



m'^'AlT'^ 



/3o^ 



'9ib{p) aibip) 



) 



N^' 



1 



(1-^)/^' 



p 



47r2 

,4 



Ar2^^'^4^2 



/^l „2 /, ^^ /'„2 ^,.^ „2 

9lbip){'^ + ^5afc(/^)) -9afe(M)(i + iralbifJ-)) 



1 p 

(^ ~ ^^)^ox;;2 (i + ^albifj')) ialbif^) - albip)) 



(h 2 



(3.32) 



The second term in the last line is a contact term, but the first term depends on gabifJ^), 
therefore it is not RGF Hence Eq.(3.32) is not exactly RGI even after subtracting the 
contact terms. The scale dependence in the physical term is due to the fact that the 
correlator is not exact but it is computed to a finite order of perturbation theory. We notice 
that the scale dependence occurs at order of g"^ only and in any case it does not afi'ect the 
structure of the universal UV behavior but only the overall coefficient in the RG estimate. 
Yet it is interesting to determine the precise overall coefficient of the asymptotic behavior. 
This is done for the correlator of l-^trF"^ in SU{3) QCD in sect. (3.7) by assuming its 
i?G-invariance, instead of checking it to a finite order of perturbation theory as we just did. 



3.3 Correlator of ^trF^ in SU{N) YM (three loops) 



We now present the result for the correlator of j^trF'^. We recall that in this case we do not 
expect to get a RGI function to all orders in perturbation theory. We start from Eq.(3.27) 
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and we write it as: 



£trF\p)^\F\-p)) 

iV iV f 



1 ^P' _4... 1 (l + |5i(A^)) 



X 



2 2 : 

1 - l^oglbifJ-) log — - 2/3i5ab(/^) log — + PhtbifJ-) log^ - 



/"^ 



M^ 



/^ 



-(l-^)|;2 5a6(/^)^ 



Af2M7r2 



T) T) T) 

1 - /305ab(/^) log -^ - /3l5afe(M) log -2 + Phabil^) log^ -f 

A* A* /^ _ 



(3.33) 



Surprisingly we notice that the term that depends on the product gabip)gab{p), that is not 
RGI, is of the same order of g'^ as the non-RGI terms in the correlator in Eq.(3.32), that 
must be RGI. 

3.4 Correlator of ^trFF in SU{N) YM (three loops) 

2 ~ 
We repeat the same steps to find the i?G-improved expression for the correlator of j^trFF, 

that is RGF The three- loop correlator reads: 



2 2 

(—trFF(p)—trFF(-p)) 



l-l^)7r2 5M5(/^)log-2 



JV2 ) 47^2 ^MS 
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Now we perform a generic change of scheme as in Eq.(3.23): 

alb = 5|?5(/^)(i + "^If5(^) + ^5fe(^)) 
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(3.35) 
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The correlator becomes: 
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Again we impose the conditions to eliminate the finite terms: 
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With this choice of a the coefficient of the g'^ log ^ term becomes: 



/2 + 3/3oa = /2 + -/o/3o 
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Substituting in the correlator we get: 
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We notice that the expression in square brackets is the two-loop Z factor implied by the 
anomalous dimension of trFF computed in Eq.(2.39). It renormalizes two powers of g{p). 
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Therefore, the correlator reads: 
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The second term in the last line is scale dependent to the order of 5^ as the term that 
occurs in the correlator of j^trF"^ in Eq.(3.33). 

3.5 Correlator of ^trF'"^ in SU{N) YM (three loops) 

We now sum the two results for the correlators of trF"^ and of trFF to obtain the correlator 
of trF~ . Indeed, we recall that: 

{trF-\p)tTF-\-p))^^^^ = 2 {tTF\p)tTF\-p))^^^^ + 2 {trFF{p)trFF{-p))^^^ (3.41) 

Summing the two results in Eq.(3.40) and in Eq.(3.33) we obtain: 

£trF-\p)^\F"\-p)) = 

iV iv conn 

= (1 - m^h^J^^^^-bip) + 91,{P) - 9lb{l^) - 9i{p) + ^9afe(/^) 
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-|5i(P) + 0(/)) (3.42) 

where surprisingly the mixed terms 5'^(p)fi'^(^) cancel to the order of g^ . There is no 
perturbative explanation for such cancellation, but conjecturally the cancellation occurs 
because of the RG invariance of the non-perturbative formula Eq.(1.8) in the TFT for 
the L = 2 ground state [4] [5] of the large- A^ one-loop integrable sector of Ferretti-Heise- 
Zarembo (see sect. (4)). We can express the last result in terms of the coupling constant in 
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the MS scheme: 
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(3.43) 



that coincides with Eq.(1.4). 

3.6 Scalar correlators in SU{3) QCD with n; massless Dirac fermions 

In this section we derive the i?G-improved expression for the correlators of trF"^ and of 
iiSXMltrF^ in QCD. 

The three-loop perturbative result for the imaginary part of the correlator of trF'^ in 
QCD with 71; massless Dirac femions is [14]: 
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We write the correlator in terms of the coupling gvM in the MS scheme instead of Ug : 
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If we suppose the correlator to be of the form: 
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its imaginary part is: 
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Comparing Eq.(3.47) and Eq.(3.45) we get: 
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Now we repeat the same steps as in the ni = case. We change renormalization scheme in 
order to cancel the finite parts: 
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We use the perturbative expression for the renormalized coupling constant with two-loop 
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accuracy: 
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where the tilde refers to the QCD coefficients of the (3 function: 
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We consider now the correlator of gyj^trF^ 
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Choosing n = -# to cancel the finite term of order of g'^ in the square brackets we get for 
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as predicted by Eq.(2.35) and by the computational experience gained in the pure YM 
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case. To cancel the finite term of order of g^ we put: 
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Therefore, the correlator reads: 
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Now we follow the same steps as in the n; = case. The only differences are the coefficients 
of the /3 function and the parameters u, v that define the new renormalization scheme. The 
result is: 
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We recall that u = -^, therefore: 
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Now we multiply the RHS of Eq.(3.58) by °^^^ ^ . Indeed, this is necessary to take into 

\ guv J 
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account the change of scheme performed to compute Eq.(3.58). The additional factor is: 
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Therefore, the correlator in Eq.(3.58) becomes: 
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that has some dependence on the scale fi even after subtracting the contact terms. In the 
next section we get rid of this dependence by using a different method, that assumes the 
RG invariance of the correlator instead of checking it. 

In any case the universal UV asymptotic behavior is in agreement with the RG esti- 
mate, i.e.: 
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3.7 i?G-invariant scalar correlator in SU{?>) QGD with ni massless Dirac fermions 

Firstly, we check the correctness of the finite parts of the scalar correlator in QGD, recon- 
structed in sect. (3.6) from its imaginary part, thanks to another result reported in [30]: 
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with: 



/n(p) = {trF'^{p)trF^i-p)) 



(3.63) 



where we have changed the overall normalization factor of the correlator with respect to 
[30] to be coherent with the one used in this paper. We perform the derivative at p'^ = fj? 
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of the correlator obtained in sect. (3.6): 
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We recall that: 
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It is easy to verify that Eq.(3.62) and Eq.(3.65) are in agreement. Indeed: 
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From Eq.(3.64) it follows the derivative of the correlator of B^SXMli^p'^ with two-loop 
accuracy: 
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Secondly, we write gYAiip) instead of gYA^ifJ-) in Eq.(3.67) since log ^ = ^ p"^ = fj? 
in order to get the large-momentum correlator in a manifestly RGI form. Exploiting the 
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definition of the /3 function we can express dlogp^ in terms of dg{p): 
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Eq.(3.69) gives the manifestly RGI form of the correlator after subtracting the /^-dependent 
contact terms. 

3.8 Correlators in the coordinate representation 

In this section we find the i?G-improved expression for the perturbative correlators in the 
coordinate representation. This procedure has the main advantage that in the coordinate 
representation the contact terms do not occur, since they are eliminated by the Fourier 
transform. Indeed, the Fourier transform of p'^ is: 
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that is supported only at x = 0. This implies that at points different from zero the contact 
terms do not occur. The RG improvement and the Fourier transform must commute up to 
perhaps finite scheme-dependent terms. Therefore, in this way we get another check of the 
asymptotic behavior. In appendix A we compute the Fourier transforms necessary to pass 
from the momentum to the coordinate representation. In particular we use the following 



34 



results: 
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Using these formulae to compute the Fourier transform of the two-loop perturbative result 
in Eq.(3.13) we get, disregarding the finite parts in Eq.(3.13): 
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Firstly, the Fourier transform has produced a new finite part. Secondly, the coefficient 
of the logarithm in the square brackets is multiplied by two after the Fourier transform. 
This implies that the factor in the square brackets renormalizes four powers of ffMsifJ'), as 
opposed to the momentum representation, where only two powers of the coupling constant 
were renormalized. This is as expected, since in the coordinate representation the correlator 
is multiplicatively renormalizable as implied by Eq.(2.41). 

To eliminate the finite term arising from the Fourier transform we change scheme 
defining: 
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where g{x) is the one-loop running coupling in the coordinate scheme [29]: 
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Therefore, the renormahzation group improved one-loop asymptotic expression for the cor- 
relator is: 

-fl-^ r„ „_ . . (3.76) 
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The Fourier transform provides automatically the change in sign necessary to obtain a 
positive expression. This is due to the fact that in the coordinate representation contact 
terms do not occur. We now go one step further performing the Fourier transform of the 
three-loop propagators in Eq.(3.8) and in Eq.(3.10). We start with the scalar correlator up 
to the overall normalization: 
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The following scheme redefinition: 
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cancels the finite term of order of g^ in the square brackets and some terms of order of 
9^ log ^5-^7 leaving only the term proportional to — 4/3i. Moreover, the finite term of order 
of 5^ in the square brackets is cancelled by a suitable choice of t, as in the previous section. 
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Eq.(3.77) now reads: 
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The scale dependent term in Eq.(3.79) occurs now at the order oi g^ , while in the momentum 
representation occurred at the order of g"^. Now we multiply both sides of Eq.(3.79) by 
(1 + ggstif^)) ) i-6- by the factor necessary to make the correlator RGI. Reinserting the 
overall normalization, we obtain: 
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As a result the possible scale dependence is of order of g^ . 

Performing the same steps for the pseudoscalar correlator in Eq.(3.10) we get: 
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The correlator in Eq.(3.81) is iZG/ in the coordinate representation with three-loop ac- 
curacy, while in the momentum representation scale-dependent terms of the order of g'^ 
occurred in Eq.(3.40). As in the momentum representation we find the correlator of trF~ 
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summing the double of the scalar Eq.(3.79) and pseudoscalar Eq.(3.81) correlators. We 
obtain: 
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The scale dependence enters the term of order of g^ as in the momentum representation in 
Eq.(3.42). 

We check the correctness of the separation of the contact terms performed in the 
momentum representation. We verify to the order of the leading logarithm that the Fourier 
transform of the i?G-improved expression in the momentum representation in Eq.(3.43) 
is equal to Eq.(3.82) in the coordinate representation. Within the leading logarithmic 
accuracy it is sufficient to put g'^{p)'. 
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Therefore, the Fourier transform of the correlator in Eq.(3.43) can be computed reducing 
it to a series of positive powers of logarithms: 
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We extract the leading logarithms of this Fourier transform. By leading we mean terms 
that have the highest power of logarithm with the power of g fixed. We use Eq.(A.20) that 
furnishes the leading logarithm of the Fourier transform: 
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We compare it with the ASD correlator in the coordinate representation Eq.(3.82): 
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We want to prove that the two series in Eq.(3.86) and in Eq.(3.87) are equal. The proof is 
by induction. We prove it for the first non trivial term, i.e. for / = 1 in Eq.(3.86): 



'- AT2) 



1.1 2^ 



iV2/vr2/3o7r2x* 



■r(4)/(M) 



(3.88) 



that is equal to the term obtained from Eq.(3.87) putting n = / = 0. Assuming that the 

m—l 



equality is valid up to the order of I log 



log- 



X II 



we show that it holds at the order of 



. Indeed, the ?7i-power of the logarithm occurs in Eq.(3.86) for I = m + 1: 



1 . 1 



P 



9\^^) 



N^'^'PoJ ' l + /3o<72(^)logg^ ("^ 



^ip-x 






f:i-iy'Vwl^9''if^)'-^^-Jiog ' 



1=0 



vr^ x° 



9 9 



l-l 



+ 



+ (-i)'"^^^^^i^/5o-r(4),-+^(^)^iog ^,^, 



(3.89) 



The ?Ti-th power of the logarithm in Eq.(3.87) occurs for the m+1 couples (n, /) such that 
I + n = m: 



(i-i^)z^ir(4)5^(-) 



iV2^vr2/3o7r2a;8 



1 . 1 



N"^' tt'^Po 



n+l<.m~l 



tt'^x^ 



r(4)/(^) E (-l)"+'(/3o5'(M)log 



4 -sn+l 



n,l=0 



x2^2^ 



+ 



n+l=m 



nJ=0 ^ 



1 . 1 



N"^' TT^pQ 



r)5 n+l<.m—l . 

n,l=0 



tt'^x^ 



x'^fl'^' 



+ (-l)™^^^^i^/5o"T(4),--^^(^)( lo. ^,,,, 



TT^X" 



x'/.r 



(3.90) 



For the inductive hypothesis the first term in the last expression is equal to the first one in 



39 



Eq.(3.89), i.e.: 



(1--) — 



n+l<m—l 2 2 



n,/=0 



1 1 r ™ 7 -n^^NnS 1 / A \ /-I 

^ iV2^7r2/3o 



E(-l)'-.^(.)/3^.^'(.)^^;^(log;;^) 



i=0 

^m 4 

by induction is complete. 



(3.91) 



The remaining terms, i.e. the terms of order of log'" -2—3-, are equal and therefore the proof 



4 ASD correlator in the Topological Field Theory 

We briefly summarize the results for the glueball propagators in the TFT underlying large- 
TV YM [2] [3] [4] [5] [6]. For the ASD glueball propagator [5] [4] ^r 

22 ..00 iL2 4 \ 6 

(|-(^-(p))f,-{^-(-p))>_ ^ ^ E j;4^ + ... (4.1) 

Besides, in the TFT the two-point correlators of certain scalar operators O2L of naive 
dimension D = 2L that are homogeneous polynomials of degree L in the ASD curvature 
-F~[5][4] can be computed asymptotically for large L: 






(02l(p)02l(-p)>,_ = constY, '^ p2%A^^ (4-2) 



The operators O2L occur as the ground state in the integrable sector of large- A^ YM of 
Ferreti-Heise-Zarembo [28] asymptotically for large L. Ferreti-Heise-Zarembo have com- 
puted their one-loop anomalous dimension for large L [28]: 

The ground state for L = 2 is the ASD operator that occurs in Eq.(4.1) for which 70(04) ~ 
2/3o exactly. 

In Eq.(4.1) and in Eq.(4.2) A^ is the RG invariant scale in the scheme in which it 

coincides with the mass gap. The functions g'^(-r^) and Z{-^) are the solutions of the 

w w 



^We use here a manifestly covariant notation as opposed to the one in the TFT [5] [4]. 



differential equations: 

dg _-/3o/+(4^5^^ 



P 



dlogp l-^g^ 

dlogZ ^ 

— = 2705 + • • • 

alogp 

where p is equal to the square root of p^. The definitions of g;^ and Z^ are: 

9k = 9{k) (4.5) 

Zk = Z{k) (4.6) 

In [2] it is shown that Eq.(4.4) reproduces the correct universal one-loop and two-loop 
coefficients of the perturbative j3 function of pure YM. Indeed, substituting in Eq.(4.4) we 
get: 



dg _ -/3o5^ + (1^9^ 



51ogP l-(4^ff2 



+ 



o 3 , ^70 5 4/3o 5 

-/3o5 + TiTToa - TiTToa +■■■ 



(47r)2^ (47r)2- 
/3o5' - /3i/ + • • • (4.7) 



-^°^' + (4^y^'-|(4^^' + 



where: 



These are the correct one- and two-loop coefficients that arise in perturbation theory of pure 
YM for the 't Hooft coupling. Therefore, the renormalization-group improved universal 
asymptotic behavior of g^ is: 



.2 1 / /3i loglogf \ ^/ 1 

/3ologfV Pi logf ) Vlog'f 



^i-TT^^ 1-^^^^FI^ +0 r4T (4.10) 
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and the renormalization group improved universal asymptotic behavior of Zf, is: 

In this section we find the asymptotics of the ASD propagator in Eqs.(4.1) and of the large- 
L propagator in Eq.(4.2) at the order of the leading and of the next-to-leading logarithms 
following the technique employed in [5] at the order of the leading logarithm. 

To find the asymptotics of the glueball propagator for large L in Eq.(4.2) we follow the 
strategy explained in sect. (1.4) for the ASD correlator. Firstly, we highlight the physical 
terms contained in Eq.(4.2) neglecting the non-physical contact terms. Secondly, we extract 
the asymptotic behavior writing the sum in Eq.(4.2) as an integral [5]. Finally, we use the 
leading and next-to-leading expression for Z^, in Eq.(4.11) to compare Eq.(4.2) with RG- 
improved perturbation theory. 

We write Eq.(4.2) as [5] [4]: 



E 



((fcA|, + p^)(fc A|,-p^) +p^)^-iZ-^A|, 
p2 + /,A2 



fe=l ^ ' w 

. 4L-4 V^ k ^^W 
00 L— 1 



fc=lm=l ^ ^ 

00 y— L A 2 

^ 4L-4y^ k w ^___ (4,12) 

where the dots stand for contact terms. 

As in sect. (1.4) we use the Euler-McLaurin formula to approximate the sum to an 
integral [5] [4]: 

E Gk{p) = / G,{p)dk -Y,^ di~'Gk{p) ^_^ (4.13) 

fc=fci -^^1 j=i ^- "-"^ 

In our case the terms proportional to the Bernoulli numbers involve negative powers of p 
and they are therefore subleading with respect to the first term, hence we ignore them. 
We obtain: 

E3iTif~[7^* ("4) 



,.,p' + kAi^ ii t + jj- 



w 
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In order to compare Eq.(4.14) to the i?G-improved perturbation theory, we substitute for 
Z^ its leading and next-to-leading logarithmic behavior given by Eq.(4.11). We define: 



and: 



/ 70 J 
Po 



A2_ 
w 



(4.15) 



(4.16) 



The integral that determines the leading asymptotic behavior is: 



lliv^ 



1 



dk 



,/3olog(f); fc + z^ 
The next-to-leading logarithmic behavior is determined by: 

' /3iloglog(f)^ 



(4.17) 



itiy) 



1 



dk 



(4.18) 



/olog(f) V ^0 log(f) jj k + u 

7' = 2 for the ASD correlator and 7' = 3^L for the large-L correlator. We show in the 
following that the leading and next-to- leading behavior of ll^f) is: 

/3i 



It{^) 



1 

IqL - /3o 



1 

/3olog^ 

w 



log log -^ 



1 



w 



10 r„i 



(4.19) 



Therefore, the asymptotic behavior of the correlator of the TFT for large-L is: 

1 , o, .^3lL-l 



4L-4 



{02L{p)02L{-p)),onn-P 



It agrees with the naive RG estimate Eq.(2.43). 



loL - /3o 



[g\p))-o' 



(4.20) 



4.1 Asymptotic series to the order of the leading logarithm 

We now perform an explicit expansion in series of ll{i')- Firstly, we change variables from 
k to k + V. 



Ili^) 



1 



i+A/3olog(V)/ ^ 



dk 



We have that: 



[log(^)]^^' = [log(^)]-^' 



1 + 



log(l - ^) 
log(l) 



It is easy to see that if c < 1: 



log(l - i 
log(l) 



< 1 



(4.21) 



(4.22) 



(4.23) 
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We define: 



log(^) 



(4.24) 



and we exploit the binomial formula [31]: 



{i+ery = Y^i^'^'-'\-iYe^ 



r=0 



(4.25) 



to obtain a series expansion. We proceed order by order in e. At the order of e^ the only 
contribution is: 

(4.26) 



7e 



e can be further expanded in powers of r/ = p-, since in the integration domain 77 < 1: 



log(l - r/) = ^ !]■■ 



m=l 



m 



Up to the order of r/^ this expansion reads: 



(4.27) 



/ / 



k'logm 



Substituting in /^(i^) we get: 



-[/3olog(^)]"^ 
k c 



1 + 



'l+u 

foo 1 u 

~ / _[/5olog(-)]-^ 



log(l - f) 
log(l) 



l+v 

00 



A;' 



1 + 7- 



dk' 
dk' 



A;'log(^) 

1 k' , C^^ ^ I k' , 

-[/5olog(-)]-^dfc'+7Vy^ ^^/?o-^ [log(-)]-7 -irffc' 



From the first integral it follows the leading asymptotic behavior [5]: 



°° 1 k' , 1 

^[Po\og{-)r< dk' = -—^p- 



l+u 



log 



l + U 



-7'+l 



Since for large f. 



log 



l + u 



-1 



(logi/ 



-1 



it follows the leading asymptotic behavior of Eq.(4.2) [5]: 



{02l{p)02l{-p)), 



p 



4L-4 



IqL - /3o 



,/3olog^ 






(4.28) 



(4.29) 



(4.30) 



(4.31) 



(4.32) 
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Performing the same steps for the ASD correlator, i.e. for L = 2 and -y' = 2, we get: 
(/3olo: 



TT^n 



f)-^ 



k + iy 



7r2/3o 



/3olog 



1 + 5 



vr2/3o/3olog-S 



(4.33) 



w 



that agrees with the leading logarithm of the asymptotic behavior in Eq.(1.3). 

We now compute the second term in the last line of Eq.(4.29), that is the first subleading 
term. We write it as: 

-^J^^ / [log{k)r^'-'dk (4.34) 

c 

and we integrate by parts: 



,,. ^.[logikr^-'dk 



i^Po 



7' 



c 



[log(A:)] 



-y-i 



1+1^ 



r-OO JJ, 



:. fc2 






(7' + l)y^^[log(A:)]' 



Y-2 



(4.35) 



We notice that the second term in the last line has the same structure as the original integral 
but with a more negative power of the logarithm. This implies that it is a less relevant term. 
Furthermore, since performing integration by parts repeatedly we always obtain integrals 
with the same structure, we can derive a possibly asymptotic series expansion for Eq.(4.34): 



7V/3o 



7 /-oo 



dk 



[log(A;)] 



^y-i 



oo 



i + y. 



— 7 — i — s 



^ s=0 \t=0 / 






s=0 
2 oo 



jrEi-^riW+t) 



W s=0 



\t=0 



log 



w 



(4.36) 



Now that we have understood the technique, we derive a complete expression taking 
into account all the terms coming from the expansion of the logarithm in Eq.(4.27), simply 
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l-i 



substituting it in I^iy): 

ri^[/5oiog(^)]-' 

Jl+v f^ C 

= r f [/3oiog(^)] 



/ 2^ m=l mfc™ 



,A:'„_y 



v/5o-' E 



»n=l 



m 



1+v 



J^/m+1 



./ n-7'-i 



log(-) 
c 



(4.37) 



Focusing on the second term: 



7 



'/So-^' E 



771 



,/ n-7'-i 



/^o-^'y E 



(-l)z." 



m=l 



TTIC" 



log(A;)]-^'"^ 



log(-) 
c 



oo 



*"'■ E 



(-1)1-' 



m=l 

oo oo 



mc" 






1+,. ^' Jl+}± mk^+^ 



E \-'-) ^ ill 
m*+i (1 



.5 = 



/3o"^'EE(-irMT 



+ z> 






m' 



i+2 



m=l s=0 

Therefore, at the first order in e we get: 

log(l-l^)^-^' 



log 
log 



1 + u 
c 

l + u 



-7'-l-s' 



-7'-s-l 



(4.38) 



°° 1 y , 

dfc'-[/3olog(-)]-^ 



i+i^ 



/fc' 



1 + 



log(l) 



~ r^[/3oiog(-)]-^' 



l+I^ 



oo oo 



+/3o-^'EE(-i)Mt 



m=l s=0 



+ Z> 



rc=o(7'+t) 



rri' 



i+2 



log 



l + l' 



-7'-s-l 



(4.39) 



We find the subleading behavior keeping only the terms with s = in Eq.(4.39). We obtain 
in the large v limit: 



/•oo 1 y 

/ dk'-[P,\og{-) 



1 



log(l - f) 
log(l) 



- r ^[/^olog(|)]-' +7'/3o-^' [log-]-^'-^ E ^2 

•^1+'^ m=0 

~ ^/3o"^' [logz.]-^'+i + 7'/3o-T'c(2)[logi.]-^'-i (4.40) 

2 

It is interesting to notice that the transcendental function (^(2) = ^ occurs, as it often does 
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in Feynman-graph computations. 

4.2 Asymptotic series to the order of the next-to- leading logarithm 

We now perform a series expansion of /^ (^)- 



I'M 



/?o 



1 (^ /3iloglog(f)^^ dk 



/3o-" 



/3o 



^log(f) 
^ 1 



/3o^ log(f) )) k + u 
p, loglog{!^)\y' dk 



+. \log(^) \ P^ log(^) 



k 



/3o 



l+u 



1+u 






log( 
log( 

fOO 

1+v 



,k — V ^ 
c ' 

k — V , 



-^ ( ^f3i loglog{!^) \ dk 

^ Q-2 l„„/fc- 



'^ dk 

— + 



log( ) 



-y-i 



/3o^ log(^) / k 



k — u dk 
loglog(^— ) — 



(4.41) 



The first integral has been evaluated in the previous section and the second term is the new 
contribution. We evaluate it at the leading order by changing variables and integrating by 
parts: 



' 02 Pa 

Po Jl+v 



log( 

c 



n-y-i 



k — u dk 
loglog(^^) — 



7^/5o 

Po ^1+;^ 



log(-) 



-y-i 



loglog(-) — 



7'§/5o"'' /" t-^'-^\og{t)dt 
/3o ^g^ 



/ A «-y 



- (^log(^)J loglog(^) + ^ (^log(^. 



(4.42) 
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The second term in brackets is subleading with respect to the first one. Putting together 
Eq.(4.42) and Eq.(4.30) we get for ll{v) : 



/3o"' 



-/ 1 / /3iloglog(|)\V 



1 Uog(!) V /3o' log(f) 



dk 
k + v 



T^^Po' lo: 



l + Z^A^^'+l /3l.-7 



7' 



/3o~^ A 1 + ^ 



y-1 
1 



log- 



/3olo^ 



-7'+l 



1 + Z^ 



/3i(7'-l) 



c 



^P^^ (log(— )) loglog(^) 



Pi 



log I 



l + v. 



-1 



log log ( 



c 



-7'+l 



/3i 



/3o(y - 1 



c 



1 - ;^ ( i°g(^^) ) iogiog(^^) 



c 



7'-l 



/5o(7' 



1 \7' 
A- 



(4.43) 



This result agrees with the RGI perturbative estimate in Eq.(2.43). Repeating the same 
steps for the ASD correlator we get: 



TT^ Ji+u k' c 7r^/3^ 71+^ 



log( 



.k — V , 



,k — u ^dk 
loglog(^^) — 



1 



vr^A 



-g\p) + 



1 






(4.44) 



Again this result agrees with the universal behavior of the i?G-improved perturbation theory 
inEq.(3.21). 
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4.3 Link with the Lerch transcendent and the polylogarithmic function 

We may obtain the asymptotic behavior by a different method as an independent check, 
relating the relevant integrals to special functions and employing the known asymptotic 
behavior of the special functions. 

We briefly recall the definition of the Lerch Zeta function [32, 33]: 

°o 2-Ki\n 

L{X,s,a) = ^-—— (4.45) 

Setting z = e^^^\ we obtain the Lerch transcendent [32, 33]: 

oo ^ 

$(z,s,a) = ^-^— (4.46) 

The Lerch transcendent admits the integral representation: 

which is valid for Re(a) > A Re(s) > A jz] < 1 or Re(a) > A Re(s) > 1 A \z\ = 1. 
The Lerch transcendent can be analytically continued to the region [34]: 

M = {{z, s, a) G (C \ {0}) X C X (C \ Z)} (4.48) 

Moreover, we exploit the following recursive formula: 

'^{z,s,a) = z^<^{z,s,a + l) + Y,-r^^r^ (4-49) 

fc=o ^" ' 

Finally, we use the relationship between the Lerch transcendent and the polylogarithmic 
function [33, 35]: 

U,{z) = z^{z,s,l) (4.50) 

where the polylogarithmic function is defined by: 

oo ^ 
k=l 

4.4 Asymptotic behavior and polylogarithmic function 

We start performing the change of variables t = log - in the integral in Eq.(4.17): 

,1, . r [^olog(|)r2' „ .-y r ^~^'' .. ur,o^ 



Setting c = e ^ in the limit e — ?• we get the upper bound: 



Ili^) = P^'' / , T-e7^^* ^ P^'' / T-^7^^* = A'-(^) (4.53) 



logil + ^e-* -'" ./. l + 7^e- 



e" 



but the upper bound is in fact asymptotic since varying c is equivalent to a change of 
scheme. Therefore, we take the hmit e — )• in order to express l\ in terms of the integral 
representation of the Lerch transcendent in Eq.(4.47). We get: 

Il{y) = V'r(-7' + l)^{-v, -i + 1,0) (4.54) 

We now exploit the relation in Eq.(4.49) with ?i = 1, a = 0, z = —v and s = —7' + 1: 

$(-!/, -7' + 1,0)= z^{-u, -7' + 1, 1) (4.55) 

Finally, we find the relation with the poly logarithmic function: 

Iliiy) = /3o"^'r(-7' + l)Li_y+i(-z/) (4.56) 

Now we use the following asymptotic expansion of Li^ [35] : 

Li,(.) ^ g,-l„l - 2-.)(2.)=.||Mr!i ,4,57) 



to find an asymptotic expansion for Ii{iy): 



-7'+l-2i 



Il(„) = /3o"^'r(-7' + 1) y (-1)^(1 - 2i-2j)(27r)2^'^i^ 
IV ; Po V / -r J /_,\ J \ j\ J (2j)!r(-7' + 2-2j) 

We get the leading behavior of /J (z^) from the j = term in Eq.(4.58): 



Keeping also the j = 1 term we obtain: 



(4.58) 



^.'M-vT,-y..,^;i^^M:^ 



A'M~ ''°,l°yl'r +^''^°"tl'°^-l"'" <«»' 
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in perfect agreement with Eq.(4.40) since ((2) = ^. Reinserting the momentum p in the 
definition of i' the asymptotic resuh is: 



'' 14, 



f— + ^oL-[/3olog(^)]-/^^-^ 



(4.61) 



Using the same technique we find the next-to-leading logarithmic behavior of I^. Indeed, 
also in this case we obtain an upper bound putting c = e~^ and taking the limit e — t- 0: 



4'(^) = c/3o 



</3c 



-rra,. 



PI t 

^logt^^^' 

Pi t 



dt 






c + ve 
dt 



l + 7^e 



i^c-t 



It.{-) 



(4.62) 
(4.63) 



but the upper bound is in fact asymptotic since varying c is equivalent to a change of 
scheme. We now expand l1_^{y): 



Il-ei^) ~ /3, 



-r r ]_(. /SiVlogA dt 
ty' ' 



P^ t J l + ^e 



L^^-t 



(4.64) 



The first term is equal to Ii_^{i'), while the second one is the new contribution. This new 
term can be linked again to the polylogarithmic function using the relation: 



t-^'-ilogt = -^t-" 
da 



(4.65) 



a=Y+l 



We find : 



i?_,(^, -7') - ilM -i) + ^^iU{^. -«) 



Pl da 



(4.66) 



a=-f'+l 



We take the limit e — )• and we perform the derivative in the asymptotic expression of 
Ii{v., —a) in Eq.(4.58). Keeping only the leading contribution we obtain: 



9 tK 



A 



a=7'+l 



y n-i) 
° r(-7' + i) 



(log u) '^ log log V 



^-j-{[ogu)~^ log log z/ 



7' 



(4.67) 
(4.68) 
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Thus the asymptotic behavior to the next-to-leading logarithmic order is: 



Il{-) 



[/^Ologl^)]-^'^^^ o 

w /^i 



P -.-^L 



loL - /3o 



^(/3ologT^) "0 log log ^^ 



1 



IqL - I3q 



,/3olog| 



/3] 



w 
log log 



El 

A^ 
w 



Pq loe 



i^^-l 



that agrees perfectly with the RG estimate Eq.(4.43). 



(4.69) 



A Fourier Transforms 

In this appendix we compute Fourier transforms of the kind: 

for L and n positive integers. We start writing: 

2\ n on / 2 \ "^ 

^ \ r d ( P ^ 



(pY(iog^) ={^^' 



da^ \fi 



a=L 



and: 



p 



2 \ " 



1 



oo _p2 



./^v r(-Q) JO 

Substituting in Eq.(A.l) and exchanging the order of integration we get: 



{ix'^y 9" 



1 



(27r)4 da^ \T{-a) jq 
The integral on p is now Gaussian and we obtain: 

1 



OO /> „2 



a=L 



(27r)4 9a" \T{-a) jq 






t2 



(2vr)' 



:(/i^)>^ 



1 



Oa" Vr(-a) 7o 



Q = L 

2 2 

o 1 a: M 



(A.l) 



(A.2) 



(A.3) 



(A.4) 



(A.5) 
(A.6) 



a=L 



We compute the last integral reducing it to a F function by the substitution t' 



4 t 



°° o 1 .2^2 

t-'^-^e—i-^dt 



^it')-' 



-a— 3 9 9 



9 9\ — Q— 2 



r (q + 2) 



The Fourier transform now reads: 

D-x " Z' 



(^2)^+2 ^n / 



4(27r)2 aa" I F(-a) V 4 



9 9\ —a— 2 



F (a + 2) 



a=L 



(A.7) 
(A.8) 
(A.9) 



(A.IO) 



We evaluate the Fourier transform in some cases by means of Mathematica. In particular 
we are interested in the cases L = 2 with n = 1,2, 3. We obtain: 

r 2n2i p"^ in-x d^P 26.3 



.2n2 



p 



2\ 2 



vr^x^ 



.2a2 



p 



(2vr)4 
„.,. d^p 2> ■?>, 10 



(pT(log^| e--^:^ = ^(-^ + 27,,-log:;:^) 



(27r)4 7r2x8 ^ 3 



x^^jp-' 



2\ 3 



(P ) log ^ e 



/^^ 



^ip-x 



.^ d^p 2^-3. 51 
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where 7^ is Euler-Mascheroni constant. 

Moreover, again by means of Mathematica, we evaluate Eq.(A.lO) for generic L and 
n = 1 or n = 2: 
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where H[L) is the harmonic number defined by: 

and ip is the digamma function defined by: 
Inverting Eqs.(A.12-A.13) we obtain: 
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We are also interested in extracting the leading logarithms in Eq.(A.lO) in the generic 
case. We obtain the leading logarithm from the terms that contain n — 1 derivatives with 

(2 2 \ ~Oi—2 
^-^ J and one derivative of p, ^^ , since otherwise we get zero because 
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— for L a positive integer: 
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The factor of n occurs because there are n such terms performing the n-th derivative. 



The limit ( r\_"J) j 



can be easily calculated knowing that the residue of the 



a— >L 



gamma function at — L is ^j— ■ The result is: 
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Therefore, the leading logarithm of the Fourier transform is: 
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